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Nanomechanical resonators having small mass, high resonance frequency and low damping rate 
are widely employed as mass detectors. We study the performances of such a detector when the 
resonator is driven into a region of nonlinear oscillations. We predict theoretically that in this 
region the system acts as a phase-sensitive mechanical amplifier. This behavior can be exploited 
to achieve noise squeezing in the output signal when homodyne detection is employed for readout. 
We show that mass sensitivity of the device in this region may exceed the upper bound imposed 
by thermomechanical noise upon the sensitivity when operating in the linear region. On the other 
hand, we show that the high mass sensitivity is accompanied by a slowing down of the response of 
the system to a change in the mass. 

PACS numbers: 42.50.Dv, 05.45.-a 



I. INTRODUCTION 

Nano-electro-mechanical systems (NEMS) serve in a 
variety of applications as sensors and actuators. Recent 
studies have demonstrated ultra-sensitive mass sensors 
based on NEMS 

Such sensors promise a broad range of applications, from 
ultra-sensitive mass spectrometers that can be used to 
detect hazardous molecules, through biological applica- 
tions at the level of a single DNA base-pair, to the study 
of fundamental questions such as the interaction of a sin- 
gle pair of molecules. In these devices mass detection 
is achieved by monitoring the resonance frequency tuo of 
one of the modes of a nanomechanical resonator. The 
dependence of loq on the effective mass m allows for sen- 
sitive detection of additional mass being adsorbed on the 
surfaces of the resonator. In such mass detectors the 
adsorbent molecules are anchored to the resonator sur- 
face either by Van der-Waals interaction, or by covalent 
bonds to linker molecules that are attached to the sur- 
face. Various analytes were used in those experiments, in- 
cluding alcohol and explosive gases, biomolecules, single 
cells, DNA molecules, and alkane chains. Currently, the 
smallest detectable mass change is 6m ~ 0.4 x 10~^^ kg 
, achieved by using a 4 ^m long silicon beam with a 
resonance frequency ujq /2Tr — 10 MHz, a quality factor 
Q of about 2, 500, and total mass m ~ 5 x 10"-'^^ kg. In 
a recent experiment Ilic et al. lOj succeeded to mea- 
sure a single DNA molecule of about 1, 600 base pairs, 
which corresponds to Sm ~ 1.6 x 10"^"'^ kg, by using a 
silicon nitride cantilever, and employing an optical de- 
tection scheme. 

In general, any detection scheme employed for moni- 
toring the mass can be characterized by two important 
figures of merit. The first is the minimum detectable 
change in mass Sm. This parameter is determined by 
the responsivity, which is defined as the derivative of the 
average output signal {X (<)) of the detector with respect 



to the mass m, the noise level, which is usually character- 
ized by the spectral density of X (t) , and by the averaging 
time r employed for measuring the output signal X {t) . 
The second figure of merit is the ring-down time iRD, 
which is a measure of the time width of the step in X (t) 
due to a sudden change in m. 

A number of factors affect the minimum detectable 
mass 6m and the ring-down time ^rd of mass detectors, 
based on nanomechanical resonators. Recent studies 0, 
[Tsl l have shown that if measurement noise is dominated 
by thermomechanical fluctuations the following hold 

TO xQiLioT Uo ) 

where fcsT is the thermal energy, Uq is the energy 
stored in the resonator, and r is the measurement av- 
eraging time, and the ring-down time is given by 



tRD = — . (2) 

Wo 

Eq. indicates that nanomechanical resonators hav- 
ing small TO and high may allow high mass sensitivity 
(small (5to). Further enhancement in the sensitivity can 
be achieved by increasing Q however, this will be ac- 
companied by an undesirable increase in the ring-down 
time, namely, slowing down the response of the system 
to changes in to. Moreover, Eq. J^l apparently suggests 
that unlimited reduction in (5m can be achieved by in- 
creasing \Jq by means of increasing the drive amplitude. 
Note however that Eq. (^), which was derived by as- 
suming the case of linear response, is not applicable in 
the nonlinear region. Thus, in order to characterize the 
performances of the system when nonlinear oscillations 
are excited by an intense drive, one has to generalize 
the analysis by taking nonlinearity into account. From a 
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more general point of view, such a generalization is inter- 
esting because it provides some insight onto the question 
of what is the range of applicability of the fluctuation- 
dissipation theorem for nonlinear systems 16]. 

In the present paper we generalize Eqs. and Q 
and extend their range of applicability by taking into ac- 
count nonlinearity in the response of the resonator to low- 
est order. Practically, characterizing the performances of 
nanomechanical mass detectors in the nonlinear region 
is important since in many cases, when a displacement 
detector with a sufficiently high sensitivity is not avail- 
able, the oscillations of the system in the linear regime 
cannot be monitored, and consequently operation is pos- 
sible only in the region of nonlinear oscillations. Another 
possibility for exploiting nonlinearity for enhancing mass 
sensitivity was recently studied theoretically by Cleland 
[rHI |. who has considered the case where the mechanical 
resonator is excited parametrically. 

When nonlinearity is taken into account to lowest order 
the resonator's dynamics can be described by the Duffing 
equation of motion [T^ . A Duffing resonator may exhibit 
bistability when driven by an external periodic force with 
amplitude p exceeding some critical value pc- Figured 
below shows the calculated response vs. drive frequency 
LOp of a Duffing resonator excited by a driving force with 
(b) sub-critical p = Pc/'^ (c) critical p = Pc, and (d) over- 
critical p = 2pc amplitude. The range of bistability in 
the (u!n, p) plane is seen in Fig. Q] (a). As was shown in 
Ref. |l8i | , high responsivity can be achieved when driving 
the resonator close to the edge of the bistability region 
El El m where the slope of the response vs. 
frequency curve approaches infinity. Note however that 
in the same region of operation an undesirable slowing 
down occurs, namely i^D can become much longer than 
its value in the linear region, which is given by Eq. ij^J. 

The detector's performances depend in general on the 
detection scheme which is being employed. Here we con- 
sider the case of a homodyne detection scheme fT3| , where 
the output signal of a displacement detector monitoring 
the mechanical motion of the resonator is mixed with a 
local oscillator at the frequency of the driving force and 
with an adjustable phase (t>hO- In the nonlinear regime 
of operation the device acts as a phase sensitive inter- 
modulation amplifier [2^ . Consequently, noise squeezing 
occurs in this regime, as was recently demonstrated ex- 
perimentally in Ref. |25j | , namely, the spectral density of 
the output signal at the IF port of the mixer depends on 
(f'hO periodically p^ . 

To optimize the operation of the system in the nonlin- 
ear region it is important to understand the role played 
by damping. In this region, in addition to linear damp- 
ing, also nonlinear damping may affect the device perfor- 
mances. Our theoretical analysis [3 shows that insta- 
bility in a Duffing resonator is accessible only when the 
nonlinear damping is sufficiently small. Moreover, a fit 
between theory and experimental results allows extract- 
ing the nonlinear damping rate. By employing such a 
fit it was found in Ref. 12711 that nonlinear damping can 
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FIG. 1: (Color online) Response of a driven Duffing resonator. 
Panel (a) shows the bistable region in the {LUp,p) plane. The 
response vs. frequency is shown in panel (b), (c) and (d) 
for sub-critical, critical, and over-critical driving force respec- 
tively. 



play a significant role in the dynamics in the nonlinear 
region, and thus we take it into account in our analysis. 

The paper is organized as follows. In section II the 
Hamiltonian of the driven Duffing resonator is intro- 
duced. The equations of motion of the system are derived 
in section III and linearized in section IV. The basins 
of attraction of the system are presented in section V. 
The ring-down time is estimated in section VI, whereas 
the case of homodyne detection is discussed in section 
VII. The calculation of the spectral density of the output 
signal of the homodyne detector, which is presented in 
section VIII, allows to calculate the minimum detectable 
mass in section IX. We conclude by comparing our find- 
ings with the linear case in section X. 



II. HAMILTONIAN 

Consider a nonlinear mechanical resonator of mass to, 
resonance frequency wq, damping rate 7, nonlinear Kerr 
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constant K, and nonlinear damping rate 73. The res- 
onator is driven by harmonic force at frequency LOp. The 
complex amplitude of the force / is written as 

/ = -2imtjp2:np^/^e''^'' , (3) 

where p is positive real, 4>p is real, and xq is give by 



(4) 



The Hamiltonian of the system is given by 

H = Hi+ Ha, + Ha, + H,2 + H,3 , (5) 

where Hi is the Hamiltonian for the driven nonlinear 
resonator 

Hi = hcuoA^'A + ^KA^'A^'AA 



(6) 



The resonator's creation and annihilation operators 
satisfy the following commutation relation 



[A, ^t] = AA'^ - At A = 1 



(7) 



The Hamiltonians Ha2 and iJ^a associated with both 
baths are given by 



Ha2 = j dujfujja\ (a>) 02 (w) , 



(8) 



III. EQUATIONS OF MOTION 

We now generate the Heisenberg equations of motion 
according to 



= [0.H1, 



(14) 



where O is an operator and H is the total Hamiltonian 



i'^^ujoA + KA^AA + 
at 



(15) 



iiua2{Lj)-iJ-e-"'"'A , (16) 



da2 (w) 
dt 



— — — = -iuja2,{uj)-id—e I'^AA. (17) 

Using the standard method of Gardiner and CoUett 
, and employing a transformation to a reference frame 
rotating at angular frequency ujp 



A = Ce-*"^"* , (18) 
yield the following equation for the operator C 



H 



a3 



dcuhujal (to) aa (to) 



(9) 



The Hamiltonian Hc2 linearly couples the bath modes 
02 {u) to the resonator mode A 



Hc2 = ^\l - I dtJ 



(10) 

whereas H^i, describes two-phonon absorptive coupling 
of the resonator mode to the bath modes 03 (cj) in 
which two resonator phonons are destroyed for every bath 
phonon created 



Hr. 



e*'^MUta3(u;) + e'"!'' al{uj)AA 



(11) 

Both phase factors 02 and ^3 are real. The bath modes 
are boson modes, satisfying the usual Bose commutation 
relations 



[a„(w),a„(tj')] = . 



(12) 
(13) 



(19) 



where 



e = [7 + 2 (cjo - cjp) + [iK + 73) etc] C - pi/2e#p 



The noise term F [t) is given by 



(20) 



F = -i^e^^^a^e'"^"^ - i2^e"^'C'^ afe^''^"^ , (21) 



where 



4" (t)^^ [ dc^e-'"(*-*°)a2 (to^io) , (22) 

VZTT J 



4"(t) = ^ / dwe-^"(*-*°)a3(to,t^) . (23) 
\/2tt J 



Note that in the noiseless case, namely when F — 0, 
the equation of motion for the displacement x of the vi- 
brating mode can be written as 
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73 

37 \xo 



dx 
lit 



2K 



^e"'""* + C.C. 
m 



(24) 



IV. LINEARIZATION 

Let C = Cm + c, where Cm is a complex number for 
which 



6 (C'm: Cm) — , 



(25) 



namely, Cm is a steady state solution of Eq. (|19|l for the 
noiseless case F — 0. When the noise term F can be 
considered as small, one can find an equation of motion 
for the fluctuation around Cm by linearizing Eq. I|19(l 



dc 

— + Wc + Vc^ = F , 
dt 



(26) 



where 



de 

w = — 

dC 



and 



c=c„ 



7 + i (wo - Wp) + 2 {iK + 73) C'mCm , 

(27) 



V 



de 



9Ct 



(*X + 73)C2 



(28) 



c=c„ 



A. Mean-Field Solution 

Using the notation 

Cm = E^/^e'^^ , (29) 
where E is positive and (pm is real, Eq. reads 

[7 + i{uJo- ujp) + {iK + 73) E] i;i/2e»0™ ^ pi/2g»0p , 

(30) 

Multiplying each side by its complex conjugate yields 



(7 + 73^^)" + (wo - + KEf 



E^p . 



(31) 



Finding E by solving the cubic polynomial Eq. H31|) 
allows calculating Cm using Eq. l|Sn|) . 



Taking the derivative of Eq. H31() with respect to the 
drive frequency Wp, one finds 



dE _ 2{luo -LUp + KE)E 



du}„ 



where 



\w\'{i-e) 



v_ 

w 



(32) 



(33) 



Similarly for the drive amplitude p 



dE 



dp \W\^{l~e) 



(34) 



Note that, as will be shown below, the value C = 1 
occurs along the edge of the bistability region. 



B. The bifurcation point 

At the bifurcation point, namely at the onset of bista- 
bility, the following holds 



'dE ~ 'dFP ~ 



(35) 



Such a point occurs only if the nonlinear damping is suf- 
ficiently small [18| , namely, only when the following con- 
dition holds 



\K\ > 7373 



(36) 



At the bifurcation point the drive frequency and ampli- 
tude are given by 



{ujp - LOo) = 7 



K 

\K\ 



8 i\K^+il) 



3V3(|i^|-V373)' 
and the resonator mode amplitude is 



(37) 



(38) 



Er 



27 



V3(|i^|-V373 



(39) 



V. BASINS OF ATTRACTION 

In the bistable region Eq. (|25|l has 3 different so- 
lutions, labeled as Ci, C2 and C3, where both stable 
solutions Ci and C3 are attractors, and the unstable 
solution C2 is a saddle point. The bistable region A 
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in the plane of parameters {ujp,p) is seen in the col- 
ormap in Fig. ^ (a). The Kerr constant in this ex- 
ample is K/ujQ = 0.001, and the damping constants are 
7/a;o = 0.02, 73 = O.IK/VS- The color in the bistable 
region A indicates the difference jCal^ — |Ci|^. The bi- 
furcation point at LOp ~ luq — [cUp — loq)^ and p = Pc is 
labeled as Ac in the figure. 

Figure |21 (a) shows some flow lines obtained by inte- 
grating Eq. (|19|l numerically for the noiseless case F = 0. 
The red and blue lines represent flow toward the attrac- 
tors at Ci and C3 respectively. The green line is the 
seperatrix, namely the boundary between the basins of 
attraction of the attractors at Ci and C3. A closer view 
of the region near Ci and C2 is given in Fig. El (b) . This 
figure shows also, an example of a random motion near 
the attractor at Ci (seen as a cyan line). The line was 
obtained by numerically integrating Eq. (|19|l with a non 
vanishing fluctuating force F. The random walk demon- 
strates noise squeezing (to be further discussed below), 
where the fluctuations obtain their largest and smallest 
values along the directions of the local principle axes (see 
appendix). 

VI. RING-DOWN TIME 

The solution of the equation of motion (|26|l was found 
inRef. [H 

/OO 
dt'G{t^t')r(t') , (40) 
-00 

where 



r (t) = + W*F {t) - VF^ (t) . (41) 

The propagator is given by 

G{t) = u{t)^ (42) 

where u(t) is the unit step function 

r 1, t>0 
u{t) = <^ 1/2, t^O , (43) 

[0, t<0 

and Ao and Ai are the eigenvalues of the homogeneous 
equation, which satisfy 
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FIG. 2: (Color online) Flow lines obtained by integrating Eq. 
H19|l for the noiseless case F = 0. The points Ci and C3 
are attractors, and C2 is a saddle point. The green line is 
the seperatrix, namely the boundary between the basins of 
attraction of both attractors. Panel (a) shows a wide view, 
whereas panel (b) shows a closer view of the region near Ci 
and C2. The cyan line near the attractor Ci in panel (b) 
demonstrates random motion in the presence of noise. 



where W is the real part of W. Thus one has 

Ao,i = T^' (1 ± ^1 + ^(C^-1)) , (46) 

or 

Xo,i=j + 2-f3E±^JiK^+-fl)E^-iujo-LUp + 2KEf . 

(47) 

We chose to characterize the ring-down time scale as 



'o + ^^=^W', (44) ^HD = (AoA,)-/^ = — ^ . (48) 

Note that in the limit C ^ 1 slowing down occurs and 
AqAi = \ W\'^ — , (45) tRD 00. This limit corresponds to the case of operating 
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the resonator near a jump point close to the edge of the The responsivity R of the detection scheme is defined 
bistabihty region. as 



VII. HOMODYNE DETECTION 

Consider the case where homodyne detection is em- 
ployed for readout. In this case the output signal of a 
displacement detector monitoring the mechanical motion 
is mixed with a local oscillator at the same frequency as 
the frequency of the pump tOp and having an adjustable 
phase (f)ho (0LO is real). The local oscillator is assumed 
to be noiseless. The output signal of the homodyne de- 
tector is proportional to 



'C(t) 



>Ct it) 



(49) 



For the stationary case of a fixed mass m the time vary- 
ing signal X^^^ (t) can be characterized by its average 



Xo = {X^,o (t)) , (50) 
and by its time auto-correlation function 

K {f - t) - {t) ~ Xo] [X^,o (^0 - ^o]) . (51) 

The correlation function is expected to be an even func- 
tion oft' — t with a maximum sA t' — t — 0. The correla- 
tion time characterizes the width of that peak. Consider 
a measurement in which X^^o (^) continuously moni- 
tored in the time interval [0,t]. Let X^ be an estimator 
of the average value of X^^^ {t) 

Xr = - r dtX^,,,{t) . (52) 

T Jo 

Clearly Xr is unbiased, and its variance is given by 



[Xr-Xof) = ^ dtj^ At' K{t'-t) . (53) 



Assuming the case where the measurement time r is 
much longer than the correlation time. For this case one 
can employ the approximation 



'(X,-Xo)') = i j dtK{t) , (54) 
or in terms of the spectral density -P^lo (^) of X^^^ (t) 



[Xr - Xo) 



27r 



Pc^.o (0) 



(55) 



R = 



dXo 



dm 



(56) 



Using Eq. (|55|l one finds that the minimum detectable 
change in mass is given by 



6m = R-^ 



1/2 



Pi!" (0) 



(57) 



Moreover, since loq is expected to be proportional to 
TO~^/^ one has 



Sm _ 2 /27r' 
m ujo \ T ^ 



^1/2 


dXo 




dujo 



pZ (0) 



(58) 



VIII. SPECTRAL DENSITY 

To calculate the spectral density P^lo W) of X4"lo (^) 
it is convenient to introduce the Fourier transform 



c(t) = 



1 



/27r ]- 



'271 J- 



Aloc {lj) e 



AloT {lj) e 



(59) 



(60) 



Assuming the bath modes are in thermal equilibrium, 
one finds 

(P(r)) = (Ft(r))=0, (61) 

(P(r)F(r'))=(pt(r)Ft(r'))=0, (62) 

{F (r) Pt (^')) ^ + Ai) 6 (r - r') (n^,) , (63) 

(Ft(r)F(r')) = (Ao + Ai)<5(r-r')((n.„> + l) . (64) 
where 

, > 1 



(65) 



and 13 = l/knT. 

In Ref. |18i I29| we have found that the following holds 



c(w) 



r(^) 



{—iuj + Ao)(-icj + Ai) 



(66) 



7 



where 

(rH) = (rtH) = o, (67) 

(r(cj')r(t^)) = Ni {uo) 5[uj + J) , (68) 

{T^Lj')T^{uj)) =Af* {lj)S{uj + lj') , (69) 

{r^{uj')r{uj)) + {r{uj')r\u;))=M2{uj)S{u;-u;') , (70) 
and 



TVi {uj) = 2W'W*Vcoth 



Af2 = 2W' (\W + iLu\^ + coth 



(71) 
(72) 



The frequency auto-correlation function of X^lo i'' 
lated to the spectral density P^lo i^) by 

(a;')^0.o (^)> = ^0Lo C'^) -5 (^ - u;') , (73) 
thus one finds 



g2i0i.o_/^^ (w) 


(io; -f- 


- Ao)(icj 


+ Xi){—iLu + Xo){—iuj - 


l-Ai) 






g-2»0Lo_;V;* (w) 




(— io; 


+ AS)(- 


-io; + A^)(ia; + AQ)(iw H 


hAt) 






AA2(t^) 




(icj + 


A5)(icj 




-Ai) 



(74) 



or in terms of the factors W and V 



where the phase factor (po is defined in Eq. HA7|) . 
The largest value 



[^0 (0)]max = coth ■ 



(77) 



is obtained when cos (0lo ^ (^^o) = 1, and the smallest 
value 



[^0(0)], 



1 2W' , (3hujo 
coth 



(78) 



"(1 + C)'|W^P 
when cos (0lo — ^o) = — 1- 

^. Integrated Spectral Density 

The integral over all frequencies of the spectral density 
is easily calculated by employing the residue theorem 



2t:W' coth 



f3hujQ 



AoAi(Ao + Ai) 



(79) 



Using Eqs. ll^ and one finds 



{uj — iAo)(w + iAo)(w — iAi)(w + iXi) 



1 P p . ^ , l + Ccos(0LO"0o) ^.pnwo 

/ "(4r n (t^]CLUJ = 7^ COtfi . 

(80) 

Thus, the integrated spectral density peaks and deeps 
simultaneously with P^lq (0). 



IX. MINIMUM DETECTABLE MASS 

To evaluate 5m using Eq. (|58|) the responsivity factor 
BXq /dojQ has to be determined. Consider a small change 
dujQ in the resonance frequency. Let Cm be the resultant 
change in the steady state amplitude Cm (here Cm is con- 
sidered as a c- number). Using Eqs. (|^ . (|?7|l . and 
one finds 



(75) 



X 2VU' coth 

2 

1. Spectral Density at lu — Q 
At frequency uj ~ Q one finds 



p l + 2Ccos((/>LO-0o) + C' 2W^^ /?;ic^o 



(81) 



Employing a coordinate transformation to the local 
principal axes (see appendix) and using Eq. PTT|) one 
finds 



\CmW^^ (Sloo) = m\ + \V\)i + 2 {\W\ - \V\) 77] 



where 



(76) 



0C = (Pm 



K - 7r/2 , 



(82) 



(83) 
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and the phase factor is defined by Eq. H29() . The 
inverse transformation Eq. IjASfl and Eq. HA7|I yield 



W 



cos (pc 1 sm (pc 

1 + C ^ 1-C 



(^^o) , (84) 



or 



Crr 



1-e 



{6coo) . (85) 



The change in Xq is given by SXo = e*'^^°c„ 



9wo 



Cm 






Re 


W 



1 - (e"2'0c 



(86) 



Finally, using Eqs. (|58|l . (|76|l . and (|86|l . and assuming 
the case of high temperature 



(87) 



one finds 



27r knT 



— = 2 



1/2 



.9 (0LO - 00 ) 



where Qcff = cjo/VK' is the effective quality factor, the 
function g is given by 




C=0.1, ([) -0.571 




C=0.99, (!)^=0.57i 




FIG. 3: The function g. Panel (a) shows (; (^lo — 0o) for 
the case = 0.1 and (pc = O.Svr, and panel (b) for the case 
= 0.99 and (j>c = O.Svr. Panel (c) shows the minimum value 
of the function g {(fiLO — </>o) vs. (fic for different values of ^. 



[1 + 2C 



cos 



cost 



1/2 



+ 0c) - CC0S((/)- (/)c)| 



(89) 



and 



X. CONCLUSIONS 



(90) 



In view of a comparison between Eq. and Eq. 

(|88|l we refer to the case where g < 1 as the case where 
the lower bound imposed upon the minimum detectable 
mass of a linear resonator is exceeded. The function 
9 {<f>w — M is plotted in Fig. JSJ (a) for the case C = 0.1 
and (pc — O.Stt, and in Fig. |j3Jl (b) for the case C = 0.99 
and (^c = O.Stt. For both cases values of g below unity are 
obtained in some range of 0lo • Figure © (c) shows the 
minimum value of the function g (0lo ~ 4>a) vs. (j)c for 
3 different values of C. In general 0.5 < t/min < 1 for all 
values of cfic and C, whereas, the lowest value f^min = 0.5 
is obtained in the limit C — > 1- This limit corresponds to 
the case of operating close to a jump point, namely close 
to the edge of the bistability region. 



In the present paper we analyze the performances of a 
nanomechanical mass detector. Both Kerr nonlinearity 
and nonlinear damping are taken into account to low- 
est order. The lower bound imposed upon the minimum 
detectable mass due to thermomechanical noise is gener- 
alized for the present case. The lowest detectable mass 
is obtained when the resonator is driven close to a jump 
point near the edge of the bistability region. However, 
in the same region slowing-down occurs in the response 
of the detector to a change in the mass (see Eq. (f^ ). 
limiting thus the detection speed. In general, for a given 
application the operating point can be chosen to opti- 
mally balance between the different requirements on the 
sensitivity and response time. 
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APPENDIX A: PRINCIPAL AXES 

Consider an expansion of the function Q near a com- 
plex number Z 



(Al) 



e{z + z,z* + z*) = eo + Wz+Vz* + o{\z 



where 60 = Qq[Z, Z*), and W and V are given by 
Eqs. ((23 and respectively. 
The transformation 



-ie"*' ie- 



(A2) 



represents axes rotation with angle (0 is real). The 
inverse transformation is given by 



Principle axes are obtained by choosing = 0o where 

WV* 



\wv\ 



(A7) 



Thus, using the notation 



WV \ 



1/2 



\WV\J 



(A8) 



one finds that in the reference frame of the principle 
axes the following hold 



s"^- {\W\ + \V\) 



R, = te^^^{\W\-\V\) 



(A9) 
(AlO) 



and 



Using this notation one finds 



Wz + Vz* = R^^ + RrjT] 



where 



i?5 = We 
Rrj = i {We 



Ve''!' , 
- Ve"^) 



(A3) 



(A4) 



(A5) 
(A6) 



Wz + Vz* = e''"^" [{\W\ + \V\)^ + i {\W\ - \V\) rj] . 

(All) 
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